function  [w,delta_w,sol_old]=set_support(nth,MI,bank_u_t,bank_delta_u_t,bank_u_p)

for i=1:MI-1
       sid(i)=i;
end

if nth > MI
    
while MI > 2
sid=input('list [MI-1 solution numbers] for the support index. sid= ');
%sid=[1 2 3];        %list MI-1 number of solutions to form the support
[c1 c2]=size(sid);
 if abs(MI-1-c2) > 0.5
    disp('Input of the support index sid is wrong');
 else
    sid=sid
    ans=input('Is sid correct? 1-yes, 0-no= ');
    if ans==1
     break;
    end    
 end   
end %while
end %if
for i=1:MI-1
w(:,i)=bank_u_t(:,sid(i));
delta_w(:,i)=bank_delta_u_t(:,sid(i));
sol_old(:,i)=bank_u_p(:,sid(i));
end
